DT1C  FILE  COPY  AD-A147  3b/ 


RADC-TR-84-45 
In-House  Report 
March  1984 


EFFICIENT  COMPUTATION  OF  REFLECTOR 
ANTENNA  APERTURE  DISTRIBUTIONS 
AND  FAR  FIELD  PATTERNS 


Hans  Steyskal 
Robert  A.  Shore 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTIOH  UNUMITEl 


MITEI^\ 

tfr' 


f 


ROME  AIR  DEVELOPMENT  CENTER 
Air  Force  Systems  Command 
Griffiss  Air  Force  Base,  NY  13441 


84  11  05  088 


This  report  has  been  reviewed  by  the  RADC  Public  Affairs  Office  (PA)  and 
is  releasable  to  the  National  Technical  Information  Service  (NT1S) .  At  NTIS 
it  will  be  releasable  to  the  general  public,  including  foreign  nations. 

RADC-TR-84-45  has  been  reviewed  and  is  approved  for  publication. 


APPROVED: 


JOHN  K.  SCHINDLER 

Chief,  Antennas  &  RF  Components  Branch 
Electromagnetic  Sciences  Division 


ALLAN  C.  SCHELL 

Chief,  Electromagnetic  Sciences  Division 


FOR  THE  COMMANDER 


'Rap 


JOHN  A.  RITZ 
Acting  Chief,  Plans  Office 


If  your  address  has  changed  or  if  you  wish  to  be  removed  from  the  RADC 
mailing  list,  or  if  the  addressee  is  no  longer  employed  by  your  organization, 
please  notify  RADC  (EEA  )  Hanscom  AFB  MA  01731.  This  will  assist  us  in 
maintaining  a  current  mailing  list. 


Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices 
on  a  pacific  document  requires  that  it  be  returned. 


SECURITY  CLASSIFICATION  OF  THIS  RAGE  (TWimi  Pit  Bnlmnd) 

REPORT  DOCUMENTATION  PAGE 

T  HfFSHT  number  i.  a ait 

RADC-TR-84-45 _ flr{ 

A.  TITLE  (ml  tut  Mid) 


_  READ  INSTRUCTIONS 
jgrORE  COMPLETING  FORM 

EIBIENT'S  CATALOG  NUMBER 


S.  TYRE  Or  RERORT  4  RERIOO  COVEREO 


EFFICIENT  COMPUTATION  OF  REFLECTOR 
ANTENNA  APERTURE  DISTRIBUTIONS  AND 
FAR  FIELD  PATTERNS 

7.  AUTHORS 

Hans  Steyskal 
Robert  A.  Shore 

’•.^PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Rome  Air  Development  Center  (EEA) 

Hans  com  AFB 

Massachusetts  01731 _ _ 

II  CONTROLLING  OFFICE  NAME  AND  AOORESS 

Rome  Air  Development  Center  (EEA) 

Hans  com  AFB 


In-house 

S.  PERFORMING  ORG.  REPORT  NUMRER 


|  •  CONTRACT  OR  GRANT  NUMBERS 


».  PERFORMING  OROANIIATlON  NAME  AND  ADDRESS  T AiK 

*  t  ■  _  .  .  *  ARt  A  a  WORK  UNIT  NUMBa NS 

Rome  Air  Development  Center  (EEA) 

Hans  com  AFB  PE  62702F  2305J303 

Massachusetts  01731 _ 

It.  CONTROLLING  OFFICE  NAME  AND  AOORESS  U  REPORT  DATE 

Rome  Air  Development  Center  (EEA)  March  1984 _ 

Hanscom  AFB  is  number  of  pages 

Massachusetts  01731 _ 

U  MONITORING  AGENCY  NAME  ft  AOORESSfJ/  dllfarant  from  Controlling  OtUca)  '5.  SECURITY  CLASS,  (oi  thla  toport) 

Unclassified 

15a.  DECLASSIFICATION' DOWNGRADING 
SCHEDULE 

_  N/A _ 

16  DISTRIBUTION  STATEMENT  (oi  thla  Ropoti) 

Approved  for  public  release;  distribution  unlimited. 


I  It.  DISTRIBUTION  ST.  4ENT  (of  ’  ■  abatract  antarad  In  Block  20,  It  dllfarant  iron i  Roport) 


18  Supplementary  TES 

RADC  Project  Engineer:  Hans  Steyskal,  RADC/EEA 


H  KEY  WORDS  < Contlnua  on  ravaraa  atda  If  nacaaaary  and  i dan t tty  by  block  numbar) 

Reflector  antenna 
Pattern  analysis 
Beam  scanning 
Limited  scan 

Array-fed  reflector _ 

3ft-  >A>tT»tCT  (Contlnua  on  ravaraa  aida  If  nacaaaaty  and  Idanttfy  by  block  numbar) 

^  A  computationally  efficient  method  is  presented  for  calculating  the  aper¬ 
ture  distribution  and  the  far  field  pattern  of  a  parabolic  reflector  with  an  off- 
axis  array  feed.  Idealized,  uncoupled  array  elements  and  realistic,  mutually 
coupled,  circular  waveguide  elements  are  considered.  Ray  optics  are  used  to 
determine  the  scalar  aperture  field  over  a  rectangular  grid  of  aperture  sam¬ 
pling  points.  Two  alternative  methods  for  performing  the  aperture  integration 
required  for  the  pattern  evaluation  are  compared:  (a)  based  on  summation  of 
subapertures  with  constant  amplitude/linear  phase,  and  (b)  based  on  the  Fast 


_ Unclassified _ 

S E CuRiTy  CLASSIFICATION  0*  THIS  PAGeT**#"  Data  Kntarad ) 


MCUWTT  CLAmUCATIQH  or  Tm»  EAOtfWii  Ol«  lal«^Q 

20.  Abstract  -  Contd. 


Fourier  Transform  (FFT)  and  sampling  theory.  Patterns  obtained  with  this 
method  are  compared  with  independently  computed  patterns  that  use  other  ap¬ 
proaches.  The  range  of  validity  of  the  aperture  field  integration  method  is 
discussed. 


_ Unclassified _ 

SCCuKlTr  CL  AlliflCATlO*  OF  T“,e  PAGEfWhtn  Dmtm  Enffd) 


Contents 


1.  INTRODUCTION  5 

2.  THEORY  6 

2.  1  Aperture  Distribution  6 

2.2  Far  Field  Pattern  10 

3.  APERTURE  BLOCKAGE  12 

4.  ILLUSTRATIVE  EXAMPLES  OF  COMPUTED  PATTERNS  13 

5.  CONCLUSIONS  22 

REFERENCES  23 

APPENDIX  A:  FEED  ARRAY  PATTERN  25 

APPENDIX  B:  FEED-APERTURE  MAPPING  33 

APPENDIX  C:  VAUDITY  OF  THE  APERTURE  INTEGRA¬ 
TION  METHOD  41 


Illustrations 


1.  Parabolic  Reflector  With  Feed  Array  7 

2.  Ray  Path  From  Feed  to  Reflector  to  Aperture  9 

3.  Blockage  Limited  Sidelobe  Level  L.  .  .  vs  Feed/Reflector 

Diameter  Ratio  d/D  13 

4.  Pattern  Computed  for  a  Uniform  Aperture  D=200A  With  r  =  4A  14 

5.  Pattern  Computed  for  a  Uniform  Aperture  With  Aperture 

Sample  Spacing  Increased  From  r=4X  to  r=8A  16 

6.  Pattern  Computed  for  a  Uniform  Aperture  With  Aperture 

Sample  Spacing  Increased  to  t  =  12X  16 

7a.  Pattern  in  the  u-plane  Through  the  Main  Beam  at 

u  =0.  05,  v  =0  19 

o  o 

7b.  Pattern  in  the  v-plane  Through  the  Main  Beam  at 

u  =0.  05,  v  =0  19 

o  o 

Al.  Feed  Coordinate  System  26 

A2.  Waveguide  Element  With  Aperture  Matching  Section  29 

Bl.  Ray  Trajectory  From  Feed  to  Aperture  34 

Cl.  Comparison  of  Patterns  Obtained  by  Scalar  Aperture 
Field  Integration  (Solid  Line)  and  by  Vector  Current 

Integration  (Circles)  42 

C2.  Geometry  of  Parabolic  Reflector  With  Line  Source  Feed  43 

C3.  Reflector  Edge  Geometry  for  GTD  Analysis  45 

Tables 

1.  Sidelobe  Positions  and  Amplitudes  -  Comparison  of 

Scalar  and  Vector  Field  Approaches  15 

2.  Sidelobe  Positions  and  Amplitudes  -  Dependence  on 

Aperture  Sample  Spacing  r  18 

3a.  Sidelobe  Positions  and  Amplitudes  in  u-plane  Cut 
Through  Main  Beam  -  Comparison  of  Scalar  and 

Vector  Field  Approaches  20 

3b.  Sidelobe  Positions  and  Amplitudes  in  v-plane  Cut 
Through  Main  Beam  -  Comparison  of  Scalar  and 

Vector  Field  Approaches  21 

4.  Computer  Time  for  One  Pattern  Cut  Over  20  Beamwidths  22 


Efficient  Computation  of  Reflector  Antenna 
Aperture  Distributions  and  Far  Field  Patterns 


1.  INTRODUCTION 

The  classic  subject  of  reflector  antennas  has  received  considerable  interest 
in  the  last  few  years.  The  reason  is  that  an  antenna  system  composed  of  a 
reflector  and  a  feed-array  combines  two  very  attractive  features,  namely  the 
inexpensive  high  gain  of  a  reflector  and  the  flexibility  of  an  array  for  beam 
shaping  or  scanning.  A  recent  effort  at  the  RADC/Electromagnetlc  Sciences 
Division  has  been  directed  at  developing  computational  capabilities  for  such 
antenna  systems  and  studying  a  limited  scan  technique  using  a  paraboloidal 
reflector  with  a  small  planar  array  feed. 

A  basic  component  of  such  an  antenna  study  is  a  computational  method  to  ob¬ 
tain  the  far  field  pattern  corresponding  to  the  particular  antenna  configuration  of 
interest.  Since  there  are  many  design  parameters  including  the  feed-array  loca¬ 
tion  and  orientation,  the  array  lattice,  and  the  element  number  and  excitation, 
this  method  must  be  computationally  efficient  in  order  to  avoid  prohibitive  com¬ 
puter  costs.  But,  the  relative  merits  of  a  particular  design  can  usually  be  judged 
from  a  limited  pattern  sector,  covering  the  main  beam  and  close-in  sidelobes  of 
the  copolarized  antenna  pattern.  This  suggests  an  approach  based  on  aperture 
integration,  which  is  highly  accurate  for  toe  main  beam  region;  and  on  toe  use  of 
scalar  theory,  which  is  a  significant  computational  sixhpliflcation. 
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The  subject  of  this  report  Is  the  presentation  of  this  computation  method.  We 
believe  the  method  to  have  two  novel  features,  not  previously  reported  in  die  vast 
literature  on  reflector  antennas.  These  are  the  use  of  analytic  expressions 
based  on  geometrical  optics  to  map  the  radiation  from  the  off-axis  feed  onto  the 
reflector  aperture,  and  the  use  of  a  rectangular  grid  of  aperture  sampling 
points,  which  simplifies  the  pattern  computation.  An  attractive  feature  of  this 
approach  is  that  in  the  process  we  determine  the  aperture  distribution,  which 
often  is  of  diagnostic  value.  Two  different  methods  of  aperture  integration  are 
pursued,  one  based  on  summation  of  subapertures  with  constant  amplitude/linear 
phase  excitation  and  the  other  based  on  the  Fast  Fourier  Transform  (FFT)  and 
sampling  theory. 1 

2.  THEORY 

2.1  Aperture  Distribution 

Consider  a  rotationally  symmetric  parabolic  reflector  fed  by  a  small,  planar 
array  with  arbitrary  location  and  orientation,  as  shown  in  Figure  1.  The  reflect¬ 
or  vertex  is  located  at  the  origin  of  the  x,y,  z  -coordinate  system  and  the  reflector 
axis  coincides  with  the  z-axis.  The  direction  to  the  observation  point  is  (0,  $), 

where  6  and  <t>  are  standard  spherical  coordinates.  The  reflector  surface  is 
given  by  the  equation, 

4  Fz  *  x2  +  y2,  z<D2/16F  .  (1) 

where  F  is  the  focal  length  and  D  is  the  diameter  of  the  reflector. 

The  feed  array  has  its  own  local  x'y'z' -coordinate  system,  whose  origin  is 
located  at  (x^,  yf,  zf)  and  whose  z'-axis  forms  the  angle  0z  with  the  positive 
z-axis. 

We  treat  the  antenna  in  the  transmit  mode  and  apply  to  the  feed  the  input 
power  Pj,  which  is  divided  into  radiated  power  Prad  and  reflected  power  Prep 
From  the  feed  we  trace  a  ray  tube  of  solid  angle  cross-section  dO  through  its 
reflection  point  on  to  the  aperture,  where  it  intercepts  an  area  dA.  Thus,  the 
complex  aperture  distribution  is  given  by 


1.  Bucci,  O.  M. ,  Franceschetti,  G. ,  and  D'Elia,  G.  (1980)  Fast  analysis  of 
large  antennas  -  A  new  computational  philosophy,  IEEE  Trans.  AntennaB 
Propag.  28:306. 


Figure  1.  Parabolic  Reflector  With  Feed  Array 
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where  the  complex  feed  pattern  f(fi)  is  normalized  such  that  |f(fi)[  represents 
the  directivity  in  the  direction  fi,  k  denotes  the  wavenumber,  and  sf  and  s&  are 
the  distances  along  the  ray  from  the  feed  to  the  reflector  and  from  the  reflector 
to  the  aperture,  respectively.  Note  that  p  is  a  "field  quantity"  and  the  intensity 

o 

|p  |  has  the  dimension  of  power  per  unit  area. 

The  antenna  pattern  is  easily  obtained  in  terms  of  the  aperture  distribution. 
The  aperture  directivity  in  the  observation  direction  (6, <t>)  is 
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where  V  ■  (u,v,  w)  *  (sin  0  cos  4,  sin  8  sin  cos  8)  is  the  unit  vector  in  the 
direction  of  observation  and  r  ■  (x_.  y„.  z  )  is  the  position  vector  of  die  a  per- 
ture  element  dA.  The  gain  pattern  of  the  antenna  system,  which  includes  spill¬ 
over  losses  but  neglects  reflected  feed  power,  is 

n  .  n  /JpJ _ dA  ... 


Similarly,  the  total  gain,  including  spillover  and  reflected  feed  power  is 
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Equations  (2)  and  (3)  form  the  basis  for  the  computer  evaluation  of  the  aperture 

distribution  p  and  the  antenna  pattern  G  . 

ap 

Returning  to  Eq.  (2)  we  note  that  the  feed  characteristics  are  described 
entirely  by  the  feed  far-field  pattern  f(fl).  This  requires  that  all  points  on  the 
reflector  indeed  be  in  the  far  field  of  the  feed  array.  However,  it  turns  out  that 
this  condition  is  not  overly  restrictive.  For  example,  a  100X  reflector  diameter 
and  50A  focal  length  allows  up  to  5A  feed  diameters. 

In  order  to  obtain  a  correct  value  for  the  absolute  gain  Gq  it  is  essential 
that  the  feed  pattern  f(I2)  be  very  accurately  known.  The  main  difficulty  lies  in 
determining  the  total  radiated  power  Pratj.  which  involves  integrating  the  feed 
pattern  over  4»  steradians.  In  Appendix  A  we  derive  the  feed  pattern  correspond¬ 
ing  to  two  different  feed-array  elements:  idealized  cos  ^-elements,  and  realistic 
circular  waveguide  elements. 

The  remaining  quantities  in  Eq.  (2),  the  mapping  function  dfi/dA  and  ray 
path  lengths  s^  and  sa<  are  purely  geometric  quantities,  which  depend  solely  on 
the  coordinates  (x^.y^.  zf)  of  the  feed  center  Ff  and  the  coordinates  (x&,  yR,  z&)  of 
the  aperture  point  Pft  (see  Figure  2).  Unfortunately,  they  are  nonlinear  functions 
of  these  coordinates  and  in  order  to  derive  explicit  expressions  for  them  we  must 
first  determine  the  coordinates  (*r>yr.  *r)  of  the  ray  reflection  point  P  on  the 
paraboloid.  We  use  the  fact  that  the  ray  path  Pf-Pr-Pa  obeys  Fermat's  principle 
and  determine  Pr  by  numerical  minimization.  Thus,  we  let  the  computer  search 
for  the  minimum  of  the  total  path  length 
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Figure  2.  Ray  Path  From  Feed 
to  Reflector  to  Aperture 


where 

222 
sf  =  [(xf-xrr  +  <yf-yr»  +  <zf-zr> 1 

1/2 

sa  =  l(xr  'xa)2  +  (*r  "  V*  +  (zr  ’ 

zr  =  <xJ  +  yJ)/4F 

za  =  D2/  16F  .  (7) 

This  proved  to  be  a  simple  and  efficient  approach  since  fast  minimization  sub¬ 
routines  for  the  computer  are  readily  available.  The  actual  derivation  of  the 
analytic  expression  for  the  mapping  dfi/dA  is  discussed  in  Appendix  B.  This 
concludes  the  derivation  of  the  complex  aperture  distribution. 

2.2  Far  Field  Pattern 

We  determine  the  far  field  pattern  G&  by  numerical  integration  of  Eq.  (3). 
For  this  purpose  we  approximate  the  circular  aperture  distribution  by  a  sum  of 
square  subapertures  with  constant  amplitude  and  linear  phase  distribution.  This 
clearly  results  in  a  rather  jagged  reflector  rim,  but  since  we  are  interested  in 


an  area  integral  this  is  of  no  consequence,  as  long  as  we  make  the  division  into 
squares  fine  enough  so  as  to  adequately  approximate  the  aperture  area.  In  most 
cases  this  represents  a  rather  mild  bound  on  the  square  size  t.  A  more  re¬ 
strictive  bound  originates  from  the  aperture  distribution  approximation.  This 
requires  t  to  be  small  enough  so  that  across  each  subaperture  the  distribution 
can  be  approximated  as 


P(*a.ya)»p(*M.yan)  e 


ill/  (x  -x  )  +  a  (y  -y  )) 
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where  (xan»yan>  denotes  the  center  of  the  n:th  subaperture,  and  the  x-  and  y- 
components  of  the  phase  gradient  at  this  point  are  and  .  These  compo¬ 
nents  are  taken  as  the  average  phase  slope  between  the  centers  of  the  two  sub¬ 
apertures  on  either  side  of  subaperture  n,  which  leads  to 

*xn  =  targp(xan+T,  yan)  -  argp<xan-T.  yan>]/2T 

tyn  =  lar«  p  (xan*  yan +  t)  ‘  ar*  p  (xan'  yan  ’  t)  1  /2  T  '  (9) 

At  the  reflector  edge,  where  these  expressions  do  not  apply,  we  average  only 
over  the  distance  r  between  the  centers  of  the  subaperture  at  the  edge  and  its 
interior  neighbor.  With  the  approximation  (8)  the  integrals  appearing  in  Eq.  (3) 
now  reduce  to 
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This  completes  the  derivation  of  the  far  field  pattern. 


We  want  to  comment  on  an  alternative  approach  to  the  preaent  pattern  com¬ 
putation.  Since  the  aperture  distribution  is  known  over  a  rectangular  grid  of 
sampling  points  it  is  tempting  to  use  a  FFT  to  evaluate  the  aperture  integral  on 
the  left-hand  side  of  Eq.  (10).  At  the  outset  of  this  study  we  felt  that  the  spatial 
variations  in  the  aperture  distribution  would  be  so  slow  that  a  constant  amplitude/ 
linear  phase  approximation  would  allow  a  much  lower  sampling  rate  than  an  FFT 
approach  and  therefore  would  be  computationally  much  more  efficient.  This  in 
fact  is  the  basic  idea  in  the  well-known  Ludwig  algorithm,  ^  which  uses  linear 
approximation  for  both  the  amplitude  and  phase  functions. 

At  a  later  stage  however,  we  became  interested  in  the  approach  of 
Franceschetti  et  al,  *  who  employ  the  FFT  to  compute  pattern  values  at  a  very 
sparse  set  of  directions,  roughly  one  value  per  beamwidth,  and  then  use  sinc- 
functions  to  interpolate  the  pattern  at  in-between  directions.  Since  the  pattern 
is  a  bandlimited  function  it  follows  from  sampling  theory  that  this  interpolation 
is  highly  accurate.  In  the  implementation  of  this  approach  we  imbedded  the 
sampled  original  circular  aperture  distribution  in  a  square  one,  with  zeros 
filled  in  over  appropriate  sections.  The  side  D'  of  the  square  was  a  parameter 
that  we  chose  slightly  larger  than  the  reflector  diameter  D,  since  in  the  FFT 
this  leads  to  more  closely  spaced  pattern  samples  and  improved  accuracy  in 
the  interpolation.  As  will  be  shown  below  with  the  aid  of  an  example,  this  method 
of  aperture  integration  is  about  20  times  faster  than  the  previous  method. 

A  final  comment  as  to  the  limitations  of  the  aperture  field  integration  method 
may  be  in  place.  Two  questions  can  be  distinguished:  (a)  over  how  large  an 
angular  sector  around  the  reflector  axis  and  (b)  over  how  many  sidelobes  around 
the  (scanned)  main  beam  can  the  method  be  expected  to  provide  reasonably  ac¬ 
curate  patterns? 

These  questions  are  discussed  in  Appendix  C.  In  regard  to  the  first  question 
it  is  concluded  that  out  to  angles  of  about  V 2F/*'  beamwidths  from  the  reflector 
axis,  the  aperture  integration  method  can  be  considered  equivalent  to  the  more 
accurate  reflector  current  integration  method.  Thus,  for  a  typical  reflector  with 
D  =  100X  and  F/D  =  0.5,  this  would  indicate  a  sector  of  ±  10  beamwidths  around  the 
axis.  However,  our  method  should  give  meaningful  results  over  a  sector  several 
times  larger  than  this,  since  we  can  expect  rather  graceful  accuracy  degradation. 

To  address  the  second  question  we  compare  our  method  with  the  Geometrical 
Theory  of  Diffraction  (GTD)  solution  for  the  field  diffracted  at  the  reflector  edge. 
The  result  is  that  the  two  methods  agree  over  a  region  of  about  (1/4)D/X  side- 
lobes  on  either  side  of  the  main  beam.  Thus,  this  number  can  be  considered  to 
be  a  bound  on  the  number  of  sidelobes  obtainable  by  the  aperture  integration 

2.  Ludwig,  A.  C.  (1968)  Computation  of  radiation  patterns  involving  numerical 
double  integration.  IEEE  Trans.  Antennas  Propag.  16:767. 
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method.  Although  we  feel  that  this  bound  may  be  a  bit  optimistic,  it  does  serve 
as  an  order  of  magnitude  indicator. 


3.  APERTURE  BLOCKAGE 


The  effects  of  aperture  blockage  are  readily  analyzed  since  this  only  requires 
setting  the  aperture  distribution  equal  to  zero  over  the  blocked  region,  before  the 
far-field  is  computed.  Presently,  provisions  are  available  in  the  computer  pro¬ 
gram  that  allow  circular  or  pie-shaped  blocked  regions. 

Frequently  however,  an  order  of  magnitude  estimate  of  the  blockage  effect 
is  more  convenient.  To  this  end  we  consider  the  feed  shadow  as  a  field  super¬ 
imposed  on  the  unperturbed  far  field  and  compare  the  radiation  intensities  in 
these  two  fields.  These  intensities  are  determined  by  the  power/gain  product  PG, 
and  therefore  the  blockage  limited  sidelobe  level 


r  '  ^blockage 

Slock  —  (p  •  G)  . 

aperture 


lp|  center 


Slock 


Ip 


average 


r)A 


ap 


center 


Ip  I 


(4)' 


average 


(12) 


where  d  and  D  denote  the  feed  and  the  reflector  diameter,  respectively,  n  is  the 
aperture  efficiency  and  lpcenter  /  Paveragel2  is  the  radiation  intensity  at  the 
aperture  center  relative  to  the  average  intensity. 

For  a  typical  example  we  choose  an  aperture  distribution  p  «  (1  -r2).  cor¬ 
responding  to  a  pattern  with  -25  dB  sidelobes.  In  this  case  rj  =  0. 75  and 
1  p center  1  paverage|2  =  3  and  therefore. 


Slock  =  4  (d/°)4 


(13) 


This  relation,  which  is  depicted  in  Figure  3,  shows  that  in  this  particular  case 
the  blockage  sidelobes  may  be  ignored  so  long  as  d/D  <  0. 15. 


Figure  3.  Blockage  Limited  Sidelobe 
Level  Lblockvs  Feed/Reflector  Diam¬ 
eter  Ratio  d/D 


4.  ILLUSTRATIVE  EXAMPLES  OF  COMPUTED  PATTERNS 

In  this  section  we  demonstrate  the  accuracy  and  the  efficiency  of  the  computa¬ 
tion  method  with  the  aid  of  a  few  examples.  Only  single  point  feeds  are  considered 
since  the  case  of  an  array  feed  simply  constitutes  a  superposition  of  such  cases. 
Patterns  shown  are  computed  with  the  piecewise  constant  amplitude/linear  phase 
aperture  distribution  approximation.  Computer  run  times  for  this  approach  and 
the  FFT  approach  are  compared  in  a  final  paragraph. 

As  a  first  case  we  consider  a  focal -point  fed  reflector.  The  feed  has  a 
l/(l+cos0)  field  pattern  which  results  in  a  uniform  reflector  aperture  distribution 
and  a  far  field  pattern  ideally  proportional  to  2J1(x)/x,  where  x  =  (kDsin  6) [2.  As¬ 
suming  a  reflector  diameter  D  =  200X  and  aperture  sample  spacing  t=4X  we  com¬ 
pute  the  pattern  shown  in  Figure  4.  As  can  be  seen,  the  sidelobes  closely  follow 
the  theoretical  pattern  envelope,  which  is  shown  by  a  dashed  line.  The  differences 
from  the  theoretical  values  are  shown  in  more  detail  in  Table  1,  column  "Scalar 
Field".  Down  to  the  -30  dB  sidelobe  level  the  differences  are  less  than  about 
0.2  dB,  they  increase  to  about  1.5  dB  at  the  -40  dB  level.  Agreement  between 
theoretical  and  computed  aperture  directivity  is  equally  good,  the  values  being 
55.96  and  55.99  dB,  respectively.  This  accuracy  was  considered  satisfactory  and 
therefore  no  patterns  were  computed  with  finer  aperture  sampling,  which  presum¬ 
ably  would  have  given  still  better  accuracy  at  the  expense  of  longer  computation 
time. 


V>V- 


•  O' 


I  , 


13 


ft  ft 


g  -30 

* 

O 

& 

-40 


-10  -9  -8  -7  -8  -5  -4  -3  -2  -1  -0  1  2  3  4  5  6  7  8  9  10 

BEAMWIDTHS  D(u-u0) 

Figure  4.  Pattern  Computed  for  a  Uniform  Aperture  D  =  200A.  With  r  «4X 


To  evaluate  the  quality  of  our  scalar  approach  to  a  vector  field  problem  we 
also  computed  the  vector  far-field  patterns  of  this  antenna  with  an  independent 

3 

computer  code  SAM*,  developed  by  Franceschetti  and  D'Elia.  The  code  assumes 
a  feed  consisting  of  a  balanced  Huygens  source  again  with  a  field  form  factor 
1/(1  + cos  0)  and  uses  reflector  current  integration.  The  sidelobes  in  the  E-  and 
H-planes  are  indicated  by  circles  in  Figure  4  and  are  also  given  numerically  in 
Table  1.  Clearly,  for  the  on-axis  case  the  polarization  dependence  is  small, 
since  the  two  patterns  differ  by  only  0.  01  dB.  From  this  data  we  infer  that  the 
SAM  code  is  highly  accurate,  so  that  it  can  serve  as  a  reference  even  for  off- 
axis  scanned  beams. 

To  evaluate  the  accuracy  degradation  with  increasing  aperture  sample  spac¬ 
ing  r,  we  recomputed  the  patterns  using  r  *  8X  and  r  *  12X,  respectively.  The 
resultant  patterns  are  shown  graphically  in  Figures  5  and  6  and  numerically  in 

♦One  might  ask  why  we  developed  an  alternative  to  the  SAM  code.  The  two  rea¬ 
sons  are  that  (a)  its  structure  is  complicated  and  we  were  unsuccessful  in  in¬ 
corporating  an  array  feed,  despite  lengthy  phone  conversations  with  G.  D'Elia 
in  Italy  and  (b)  we  felt  that  the  scalar  approach  would  provide  acceptable  data 
with  less  computer  run  time. 

3.  Franceschetti,  G. ,  and  D'Elia,  G.  (1982)  SAM  Program.  Computer  program 
for  offset  paraboloidal  reflectors.  Material  from  Reflector  Antenna  Theory, 
Computation,  and  Synthesis  course.  University  of  Southern  California, 

Los  Angeles,  Calif. .  May  1982. 
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Figure  5.  Pattern  Computed  for  a  Uniform  Aperture  With  Aperture  Sample 
Spacing  Increased  From  t  =  4X  to  t  =  8X 
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Figure  6.  Pattern  Computed  for  a  Uniform  Aperture  With  Aperture  Sample 
Spacing  Increased  to  r  =  12 X 


Table  2.  Apparently,  the  accuracy  degrades  rapidly  even  though  the  present 
uniform  aperture  distribution  constitutes  a  particularly  simple  case,  which  is 
perfectly  resolved  by  the  constant  amplitude/linear  phase  approximation. 

In  a  second  example  we  scan  the  beam  approximately  10  beamwidths  by 
moving  the  feed  off-axis  to  a  position  xf  ■  -5.  861X,  y{  =  0,  zf  ■  99.  828X,  while 
keeping  it  pointing  at  the  reflector  center.  The  feed  pattern  and  reflector  re¬ 
main  unchanged.  Using  again  an  aperture  sample  spacing  r  *  4X  we  obtain  the 
u-  and  v-plane  patterns  (through  the  main  beam  position  uQ  =  -0.  05,  vQ  =  0) 
shown  in  Figures  7a  and  7b.  As  a  reference  we  also  show  the  sidelobe  levels 
obtained  with  the  SAM  code.  No  distinction  between  the  feed  polarization  being 
parallel  or  orthogonal  to  the  scan  plane  is  made,  since  this  results  in  differ¬ 
ences  <  0.  03  dB  in  sidelobe  level.  A  numerical  comparison  of  sidelobes  ob¬ 
tained  for  the  scalar  and  vector  field  patterns  is  provided  in  Tables  3a  and  3b. 

In  the  u-plane.  Table  3a,  where  no  depolarization  occurs,  we  note  again  good 
agreement  between  the  two  approaches,  the  difference  being  £  1  dB  for  sidelobe 
levels  down  to  about  -30  dB.  In  the  v-plane.  Table  3b,  we  note  that  depolariza¬ 
tion  becomes  significant.  If  we  compare  the  pattern  of  only  the  dominant 
fl-polarized  field  component,  we  find  that  the  sidelobes  are  consistently  lower 
than  those  of  the  scalar  pattern.  This  is  due  to  the  power  lost  to  the  4 -polarized 
field  component.  When  we  add  the  powers  of  these  two  polarizations,  as  is  done 

n  o 

in  the  column  labeled  Vector  Field,  E  J  +  (E.  |  in  Table  3b,  we  find  much 

0  v 

improved  agreement  with  the  scalar  pattern.  The  differences  are  ss  0.  5  dB  down 
to  the  -30  dB  sidelobe  level. 

Finally,  we  illustrate  the  efficiency  of  the  FFT  far-field  pattern  computation 
with  one  representative  example.  We  consider  the  same  reflector,  D  ■  200X, 
F/D*  0.  5,  as  before  with  the  beam  scanned  10  beamwidths  off-axis.  The  aper¬ 
ture  sample  spacing  t  =  4X  and  for  the  FFT  we  use  64  x  64  samples,  which  makes 
the  zero-filled  aperture  diameter  D'  ~  250X.  This  choice  for  D'  results  in  com¬ 
parable  accuracy  for  the  aperture  integration  methods  based  on  the  constant 
amplitude/linear  phase  approximation  (denoted  "subapertures")  and  on  the  FFT, 
respectively.  The  computer  run  times  for  the  two  methods  are  compared  in 
Table  4,  which  lists  seconds  of  CP-time  on  a  CDC  Cyber  750  computer.  Time 
for  the  64  x  64  FFT  alone  is  about  0.  6  sec  CP-time,  which  shows  that  50  percent 
of  the  time  is  used  for  the  pattern  interpolation. 


CM 

•*. 

a> 

/< 

«  2*ffl 


■* 

in 

CM 

00 

c» 

m 

fH 

m 

H 

o 

CM 

w 

CO 

CM 

CO 

in 

in 

i 

1 

• 

• 

• 

i 

» 

1 

CM 

CO 

CO 

CO 

05 

CM 

CO 

t- 

CO 

CO 

C- 

o 

CO 

CM 

CO 

CO 

CM 

CM 

CM 

CO 

CO 

CO 

CO 

eo 

co 

CO 

c- 

00 

CO 

00 

05 

05 

CM 

CO 

in 

CO 

C- 

00 

05 

C- 

CO 

CO 

CM 

05 

CO 

CO 

o 

o 

H 

CM 

CO 

CM 

d 

eo 

1 

1 

1 

1 

1 

1 

i 

05 

CO 

in 

-<* 

00 

CM 

00 

05 

CO 

T* 

CO 

05 

H 

H 

in 

05 

CO 

H 

CM 

CM 

CM 

CO 

CO 

CO 

CO 

CO 

CO 

C- 

C- 

00 

00 

05 

05 

05 

05 

CM 

CO 

in 

CO 

C- 

GO 

05 

o 

w* 

o 

CM 

00 

05 

CO 

00 

in 

o 

o 

o 

o 

o 

o 

H 

o 

H 

i 

1 

• 

CO 

05 

o 

CO 

CO 

CM 

05 

o 

D- 

CO 

00 

H 

CO 

CO 

05 

05 

CM 

CM 

eo 

CO 

CO 

CO 

CO 

CO 

CO 

r- 

r- 

r» 

t- 

r» 

H 

CM 

co 

in' 

CO 

C- 

00 

d 

CD 

CO 

o 

CO 

c*» 

m 

H 

in 

r- 

CO 

00 

CO 

in 

05 

o 

H 

— 

CM 

CM 

CO 

CO 

CO 

CO 

CO 

-^Ncj^incor-toc: 


i 


» 


Table  3a.  Sidelobe  Positions  and  Amplitudes  in  u-plane  Cut  Through  Main 
Beam  -  Comparison  of  Scalar  and  Vector  Field  Approaches.  Main  beam 
scanned  to  u  =0.  05,  v  =  0 


Sidelobe 

No. 

Scalar  Field 

Vector  Field 

Diff 
(A  dB) 

Position 

lD(u-u0)J 

Ampl 

<-dB) 

Position 

[D(u-u0)l 

Ampl 

<-dB) 

-8 

-9.3 

31.3 

-9.3 

31.0 

-0.3 

-7 

-8.2 

28.0 

-8.3 

29.  0 

1.0 

-6 

-7.2 

26. 1 

-7.2 

26.2 

0.  1 

-5 

-6.2 

23.2 

-6.2 

23. 1 

-0.  1 

-4 

-5.  1 

19.3 

-5. 1 

19.4 

0.  1 

-3 

-4.0 

15.2 

-4.0 

15.5 

0.3 

-2 

-2.  8 

10.9 

-2.8 

11.  1 

0.2 

-1 

-1.5 

6.8 

-1.6 

7.  1 

0.3 

main  beam 

0 

0 

0 

0 

0 

1 

4.0 

29.  1 

4.0 

29.  0 

-0.  1 

2 

5.  1 

37.4 

5. 1 

36.0 

-1.4 

3 

6.  1 

37.  5 

6.  1 

35.9 

-1.6 

4 

7.  1 

35.  8 

7.  1 

38.2 

2.4 

5 

8.  1 

40.  1 

8.  1 

38.6 

-1.5 

6 

9.  1 

38.  3 

9.  1 

40.7 

2.4 

Table  4.  Computer  Time  for  One  Pattern  Cut  Over  20  Beamwldthe 


Number  of 

Time  for 

Total  Run  Time 

Pattern  Integration 
Time 

Subaperture  FFT 

FFT  Time 

Pattern 

Points 

Aperture 

Distribution 

Subaperture 

FFT 

Reduction 

Factor 

200 

11.5 

35.0 

12.8 

_ 

23.5  1.3 

18 

5.  CONCLUSIONS 

We  have  developed  a  simple  and  efficient  computational  method  for  the  eval¬ 
uation  of  reflector  antenna  patterns.  The  method  is  useful  in  situations  where  a 
large  number  of  candidate  antenna  configurations  has  to  be  analyzed,  as  is  the 
case  when  synthesizing  an  antenna  to  a  given  design  goal  or  when  scoping  a  re¬ 
flector  limited  scan  technique. 

The  three  essential  features  of  the  approach  are  (a)  the  use  of  ray  optics  to 
map  the  off-axis  feed  onto  the  reflector  aperture,  (b)  the  use  of  Fermat's  princi¬ 
ple  to  determine  the  ray  path  from  the  feed  to  a  given  point  on  the  aperture,  and 
(c)  the  use  of  a  rectangular  aperture  sampling  grid  and  FFT  to  speed  up  the 
pattern  integration.  The  pattern  and  the  absolute  gain  are  determined  for  arbi¬ 
trary  feed  locations  and  orientations,  the  only  restriction  being  that  the  reflector 
lies  entirely  in  the  Fraunhofer  region  of  the  feed.  Expressions  have  been  derived 
for  the  far  field  feed  array  pattern  corresponding  to  two  different  feed  array  ele¬ 
ments:  idealized  cos  6-elements  and  mutually  coupled,  circular  waveguide  ele¬ 
ments. 

Patterns  computed  by  the  aperture  field  integration  method  compare  well 
with  known  analytic  results  and  with  patterns  computed  independently  with  the 
SAM  code,  which  is  based  on  reflector  current  integration.  The  validity  of  the 
aperture  integration  method  is  further  discussed  in  an  appendix,  where  it  is  com¬ 
pared  with  integration  over  the  actual  reflector  surface  and  with  the  GTD  solu¬ 
tion.  A  present  limitation  lies  in  the  scalar  field  formulation,  which  neglects 
the  cross  polarization.  However,  the  latter  can  be  easily  incorporated  using 
simple  geometrical  polarization  factors,  such  as  derived  in  Reference  4. 


4.  Kauffman,  J.F.,  Croswell,  W.F.,  and  Jowers,  L.  J.  (1976)  Analysis  of  the 
radiation  patterns  of  reflector  antennas,  IEEE  Trans.  Antennas  Propag. 
24:53.  - 
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Appendix  A 


Feed  Array  Pattern 

In  this  section  we  derive  the  scalar  far-field  pattern  f  of  the  feed  array.  The 

2 

main  difficulty  is  the  proper  normalization  such  that  |f(0,  0)  |  represents  the 
directivity  in  the  direction  (0,0).*  Two  types  of  array  elements,  idealized  cos  0- 
elements  and  realistic  circular  waveguide  elements,  will  be  considered. 

We  assume  that  the  far  field  of  the  single,  isolated  feed  element  with  excita¬ 
tion  V  can  be  expressed  as 

Eo(r,  0. 0)  =  V(eofl  §  +  eQ<t)  0)  ej0  •  e‘Jkr/r  ,  (Al) 

where  eQ0  =  eo^O*0)»  eQ<j,  =  e0 are  real  functions,  0  and  0  are  unit  vectors, 

i  represents  a  constant  phase  angle,  and  e  -^/r  represents  the  spherical  wave- 

N 

front.  Consequently,  for  an  N-element  array  with  excitation  coefficients  {Vn } ^ 
the  far  field  is 

E  =  (eQ00~  +  eQ^0)  ej*  (e‘jkr/r)  E  VR  e^  '  rn  ,  (A2) 


*Note  that  in  this  appendix  r,  0,  0  stand  for  the  spherical  coordinates  in  the  local 
feed  coordinate  system,  which  in  the  main  text  carries  a  prime  sign. 


where  r  ■  (sin0cos4,  sin0sin4.  cos0)  is  the  unit  vector  in  the  observation  direc¬ 
tion  (0,4)  and  rR  is  the  position  vector  of  the  n:th  element,  see  Figure  Al. 
From  Eq.  (A2)  we  derive  a  scalar  field  E,  which  correctly  describes  the  field 
magnitude  and  phase  but  suppresses  the  polarization,  as  required  for  our  ray 
analysis,  by  setting 


-jkr  jkr  •  r_ 

E  =  eo  (e  °/r0)  L  Vn  e  " 


where 


(A3a 


(A3b 


OBSERVATION  POINT 

(r.  »,</>) 


Figure  Al.  Feed  Coordinate  System 

The  desired  scalar  far-field  pattern  f(0, 4)  is  now  defined  by  die  conditions: 

n 

|f(0, 4)1  =  directivity  in  direction  (0,4), 


(A4 


r  ■  ?n 

argf  =  arg  I  e  EVne 


which  leads  to 


U6. 4) 


or  alternatively 

1/2 


f(0, 


4vY0 

.*)«  -p - 

*rad 


e  E  V  e 
o  n 


jkr.  rn 


where  Yq  is  the  free  space  admittance  and 


prad  "  /  Yo  lEl2r2d« 

4  v 


denotes  the  total  power  radiated  by  the  feed  array. 


Al.  IDEALIZED  ARRAY  ELEMENT 


In  this  case  the  array  element  pattern  is 


3(0. 4)  ■=  | 


Vcos  e'  0  s  ir/2 


0  >  ir/2 


and  the  only  difficulty  lies  in  the  evaluation  of  the  integral. 
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ikT  *  i* 

I=  /  |eoZ  Vne  B|  dn 
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(AG 


in  the  denominator  of  Eq.  (A5a). 

Substitution  of  Eq.  (A7)  in  Eq.  (A8)  and  expansion  of  the  square  leads  to  a 
double  sum 
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where  each  term  is  given  by 
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and  the  asterisk  denotes  comples  conjugate. 

Equation  (A9b)  can  be  integrated  analytically  and  the  result  is 
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IV, 


m  =  n 


where 
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snm  =  l(xm  ‘  Xn>  +  ^  1 


(All) 


is  the  distance  between  element  m  and  element  n,  and  denotes  the  Bessel 
function  of  the  first  kind  and  order  1.  Thus,  for  cos  0-elements  the  field  pattern 
for  the  feed  array  is 
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where  the  terms  T  are  given  by  Eq.  (A10). 
mn 


A2.  CIRCULAR  WAVEGUIDE  ELEMENTS 

In  this  case  the  feed  consists  of  an  array  of  open-ended  circular,  dielectric 
loaded  waveguides  in  a  ground  plane.  Each  waveguide  is  matched  to  its  (single 
element)  radiation  admittance  and  is  assumed  to  be  driven  in  the  dominant  mode 
by  a  matched  generator. 


The  desired  pattern  function  is  now  derived  in  the  following  way.  The 

aperture  field  of  each  waveguide  is  approximated  by  the  two  orthogonal  dominant 

modes,  one  driven  and  one  parasitic.  Ideally,  the  aperture  excitation  of  each 

guide  and,  consequently,  die  radiated  field  is  directly  proportional  to  its  driving 

generator  voltage.  In  a  realistic  array  however,  the  aperture  excitation  will  be 

perturbed  by  the  mutual  coupling  between  the  elements.  We  evaluate  this  effect 

with  die  aid  of  a  separate  computer  code  based  on  Reference  A1  with  cylinder  rad- 

ius  —«o,  which  for  a  given  set  of  N  elements  driven  with  voltages  {An}j  computes 

perturbed  driven  aperture  voltages  {V' and  reflected  mode  amplitudes  {R* 

m  n  i  ||  ^  n  l 

and  corresponding  quantities  {Vn}j  and  {Rn}^  for  the  parasitic  mode,  see 
Figure  A2. 
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Figure  A2.  Waveguide  Element  With 
Aperture  Matching  Section 


The  explicit  form  for  the  y-polarized  TEll-mode  function  d(r,  $)  =  r  +  4^$  is 
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where  a  is  the  guide  radius. 

The  far  field  E  =  E  fl  +  E  radiated  from  an  electric  aperture  field 
_  A i  o$  09 

V'  *  isA1 


Al.  Steyskal,  H.  (1977)  Analysis  of  circular  waveguide  arrays  on  cylinders, 
IEEE  Trans.  Antennas  Propag.  25:610. 
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By  comparison  of  Eq.  (A14)  with  Eq.  (A3b),  we  deduce  the  scalar  element 
pattern  to  be 
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(A  15) 


where  the  superscript  (y)  indicates  that  this  pattern  corresponds  to  the  y-polar- 
ized  TE 11 -mode.  The  pattern  for  the  x -polarized  TE 11 -mode  is 


(A16) 


which  is  obtained  from  Eq.  (A15)  by  a  simple  interchange  of  sintf  and  cos$. 

For  the  proper  normalization  of  the  feed  pattern  f  we  also  must  determine 
the  total  power  radiated  by  the  array.  This  could  again  be  done  by  integrating 
the  pattern  over  4*  sr,  but  fortunately  this  time  there  is  a  simpler  method.  The 
power  incident  at  the  waveguide  apertures  Is 


YTE11  n 
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and  the  power  reflected  is 


Pr  ■  VTEU  g  <Kl2  '  lO2'  • 

where  YTE11  denotes  the  TE 11 -mode  admittance, 
radiated  power  and  therefore. 
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The  difference  must  be  the 
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The  scalar  field  pattern  now  becomes 
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where  Prad  is  given  by  Eq.  (A19)  and  the  element  pattern 


feo 

'°  US 


(x) 
o 

(y) 


is  given  by  Eqs. 


for  x -polarized  excitation 
for  y-polarized  excitation 

(A  15)  and  (A16). 
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Appendix  B 


Feed-Aperture  Mapping 


In  this  section  we  determine  the  relation  between  the  solid  angle  cross- 
section  dfi  of  a  ray  tube,  which  emanates  from  the  feed,  and  the  area  dA,  which 
this  ray  tube  intercepts  at ‘the  aperture,  see  Figure  Bl.  This  is  a  purely 
geometrical  problem  and  we  repeat  the  defining  equations  here  before  we 


elop  the  solution. 

(1)  Paraboloid  equation 

n 

+ 

CM 

X 

(2)  Focal  length 

F 

(3)  Reflector  diameter 

D 

(4)  Given  ray  origin  point 

(xf ,  yf ,  zf) 

Known  reflection  point 

(Wzr 

Given  aperture  intercept  point 

<Vya'za 

(5)  Feed-reflector  s,  = 

(V 

-xf>2  +  (yr- 

distance  1 
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(6)  Reflector-aperture 
distance 


(x  -x  )2  +  <y -yj2  + 
r  a  r  a 


2  2 
x  +  y 
r  •’r 

4F 


1  1/2 


(7)  Total  ray  path  length  8  =  sf  +  sa  =  minimum  »  (B2) 

with  respect  to  xr*  yr»  (Fermat's  principle). 

We  also  introduce  a  local,  feed-centered  coordinate  system  which  is 

translated  relative  to  the  system  x,y,  z, 

I  =  x  -  xf 

h  =  y  -  yf 

£  =  z  -  Zj.  ,  (B3) 

since  in  this  system  the  ray  tube  cross-section  df2  is  conveniently  expressed  as 

df2  =  |  sin6  de  d$|  ,  (B4) 

where  6  and  <t>  are  the  local  spherical  coordinates.* 

The  desired  relation  between  the  solid  angle  dQ  and  the  aperture  element 
dA  is  now  obtained  as  follows.  We  consider  the  coordinates  8,  <t>  that  determine 


♦Note  that  in  this  appendix  the  spherical  angles  8,  are  different  from  those  in 
Appendix  A,  since  the  ?,  n.  S -system  here  is  rotated  relative  to  the  x.y,  z-system 


the  ray  direction,  tobe  functions  of  the  coordinates  *a*ya  at  which  the  ray  inter¬ 
cepts  the  aperture. 


6 

♦-♦ka-V  •  (B5) 

The  mapping  of  an  aperture  area  element  dx&  dy&  onto  an  element  d0  d$  is 
given  by 

do  d4  ■  J  dx&  dy&  ,  (B6) 

where  the  Jacobian 


30  d<t>  30  3 <t> 

K  ^  3xa 


(B7) 


Substituting  Eq.  (B4)  in  Eq.  (B6)  and  setting  dA  =  |  dxft  dyR  |  yields  the 
desired  relation 


This  completes  the  derivation  of  dO/dA.  The  remainder  of  this  appendix  will  be 
devoted  to  the  development  of  explicit  expressions  for  the  partial  derivatives  in 
the  Jacobian. 

According  to  Fermat's  principle  Eq.  (B2)  the  reflection  point  (xr,yr,  z^)  is 
a  stationary  point  so  that 
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For  small  position  increments  dxa,  dy&  in  the  aperture  point  the  resulting 


increments  dxr>  dyr  in  the  reflection  point  will  be  such  that  Eq.  (B9)  still  holds. 
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From  Eq.  (BIO)  we  solve  for  dxr,  dyr  in  terms  of  dxg,  dya  to  obtain 


dxr  =  o  dxa  +  0  dya 


dyr  =  y  dxa  +  6  dya  . 
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The  derivation  of  explicit  expressions  for  the  second-order  partial  deriva¬ 
tives  of  8-Bf+sa  required  for  a,  y,  d  is  lengthy  but  straightforward.  The 
result  is 
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Returning  to  Eq.  (Bll),  we  note  that  the  coordinate  increments  dxr,  dyr  on  the 
reflector  surface  also  are  accompanied  by  an  increment 


dzr  =  (xrdxr  +  yrdyr)/2F  .  (B14) 

We  have  now  traced  the  aperture  point  increments  dx&,  dy&  to  corresponding 
reflection-point  increments  dxr,  dyr<  dzr  and  it  now  only  remains  for  determine 
how  these  increments  are  viewed  from  the  feed  in  terms  of  d0,  d$.  In  other 
words,  we  need  to  determine  the  unknown  coefficients  in  the  relation 

dfl  =  aj  dxr  +  a2  dyr  +  a3  dzr 

d$  =  bj  dxr  +  b2  dyr  +  b3  dzf  .  (B15) 

Once  we  have  these  coefficients  we  can  substitute  Eqs.  (Bll)  and  (B14)  in 
Eq.  (B15)  to  obtain 


dfl  =  cx  dxa  +  c2  dya 


d$  =  c3  dxa  +  c4  dya  ,  (B 

where  the  four  new  coefficients  Cj,  Cg,  Cg,  c4  now  can  be  identified  as  the  four 
partial  derivatives  required  for  the  Jacobian.  Thus,  we  have 

if-  •  Cj  =  o(»1  +  a3xr/2F)  +  T(a2  +  a3yr/2F) 
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f|-  *  c3  *  a(b1  +  b3xr/2F)  +  y(b2  +  b3yr/2F) 


|f-  -  c4  -  ^(bj  +  Vr/2F)  +  6(b2  +  b3V2F)  • 


The  coefficients  aj.ag.ag.b^.bj.bg  are  determined  as  follows.  From  the 
standard  relation  between  spherical  and  rectangular  coordinates 


cos  6 


Vi 


■5 — 5 — 

2+n2+r 


cos  4  * 


1 

t 


(B18) 


we  obtain  by  differentiation 

dfl  =  [??dC  +  r)?dn  -  t2d?]/s2  t 


d<f»  =  [-ndf  +  ?dn]/t2  .  (B19) 

where 

t  =  Vr  +  n  .  (B20) 

In  view  of  Eq.  (B3)  we  have 

d?  =  dx,  dr)  ■  dy,  d£  =  dz  (B21) 


and  when  we  substitute  Eq.  (B21)  in  Eq.  (B19)  and  identify  the  resulting  expres¬ 
sion  with  Eq.  (B15)  we  find 
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where  the  subscript  r  indicates  that  the  value  at  the  reflection  point  is  to  be  used 


Substitution  of  Eqs.  (B3)  and  (B22)  In  Eq.  (B17)  yields 
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f^-=[-P(yr-yf)  +  6(xr-xf)i/tJ 
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In  summary,  we  have  now,  via  Eqs.  (B8),  (B7),  (B23),  (B12),  and  (B13), 
expressed  the  desired  function  dU/dA  entirely  in  terms  of  the  known  coordinates 
of  the  feed  point,  reflection  point,  and  aperture  intercept  point. 


Appendix  C 


Validity  of  the  Aperture  Integration  Method 


We  want  to  briefly  comment  on  the  validity  of  the  patterns  obtained  by  the 
aperture  integration  method.  In  particular,  we  want  to  know  over  how  large  an 
angular  sector  around  the  reflector  axis  the  method  is  valid  and  how  many  side- 
lobes  about  the  (scanned)  beam  are  accurately  predicted. 

The  first  question  may  be  answered  by  comparing  our  method  of  integrating 
over  an  equivalent  planar  aperture  with  a  more  accurate  integration  over  the  actual 
curved  reflector  surface.  This  is  done  in  Reference  Cl,  where  it  is  shown  that 
the  two  methods  are  equivalent  out  to  angles  in  the  order  of  V^F/x’  times  the 
beamwidth.  For  larger  angles  there  will  be  quantitative  differences,  although  we 
would  still  expect  good  qualitative  agreement. 

One  numerical  example  illustrating  this  point  is  provided  by  the  patterns  of 
Table  3.  There,  good  agreement  between  the  aperture  field  integration  and 
current  integration  methods  is  shown  out  to  19  beamwidths  from  the  reflector 
axis,  which  is  1.4  times  further  than  the  14  beamwidths  obtained  by  the  y 2F/X- 
criterion. 

The  conservativeness  of  this  criterion  is  further  demonstrated  in  the  next 
example,  where  it  is  exceeded  by  a  factor  of  3.  We  consider  a  reflector  with 
D=  100X,  F/D=  0.  5,  with  a  balanced  Huygens  feed  located  at  x^  =  -9.  98X,  y^  *  0, 

2^  =  49.  08X,  and  with  a  l/(l+cos0)  field  pattern  pointing  at  the  reflector  center. 

Cl.  Clarke,  R.H.,  and  Brown,  J.  (1980)  Diffraction  Theory  and  Antennas, 

John  Wiley  &  Sons,  N.Y.,  p.  196-199. 


The  nominal  main-beam  position  is  uq  =  0. 174,  vq  =  0,  corresponding  to  a  scan 
angle  of  about  17.4  beamwidths.  The  u-plane  pattern  through  the  main  beam 
obtained  with  the  aperture  integration  method  is  shown  in  Figure  Cl.  The  side- 
lobe  magnitudes  and  positions  as  computed  with  the  current  integration  method 
(SAM)  are  indicated  by  circles  in  the  same  figure.  The  maximum  difference  is 
2. 7  dB  and  clearly  the  two  methods  produce  very  similar  patterns  out  to  at  least 
27  beamwidths  from  the  axis.  The  V 2F/\'  -criterion  yields  only  1C  beamwidths 
in  this  case.  The  corresponding  v-plane  patterns  through  the  main  beam  agree 
with  the  same  accuracy  as  in  the  u-plane  pattern  when  the  powers  of  the  8-  and 
d -polarized  components  are  combined.  Computations  with  two  orthogonal  feed 
polarizations  showed  that  the  patterns  are  very  insensitive  in  this  regard,  since 
the  sidelobe  peak  variations  were  less  than  0.  01  dB  down  to  the  -40  dB  sidelobe 
level. 
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Figure  Cl.  Comparison  of  Patterns  Obtained  by  Scalar  Aperture  Field 
Integration  (Solid  Line)  and  by  Vector  Current  Integration  (Circles) 


To  address  the  second  question  we  compare  the  diffraction  patterns  computed 
by  the  aperture  integration  method  and  Geometrical  Theory  of  Diffraction  (GTD), 
respectively.  The  latter  method  is  used  as  a  reference,  since  it  correctly 
describes  the  far-out  sidelobes.  Such  a  comparison  is  quite  illustrative  and  will 
be  outlined  below,  since  rather  surprisingly,  we  were  unable  to  find  it  explicitly 
in  the  literature. 

We  consider  the  2-D  case  of  a  parabolic  cylinder,  diameter  D  and  focal 
distance  F,  illuminated  by  a  focal  line  source.  The  geometry  and  notations  are 
shown  in  Figure  C2.  Note  that  the  reflector  rim  is  not  necessarily  a  knife  edge, 
but  a  general  wedge  of  angle  The  feed  has  a  symmetric  y-polarized  field 
pattern  such  that  E(s ,</)  =  E(s,  -< /).  If  we  trace  a  ray  tube  of  angular  cross- 
section  di/  from  the  feed  to  the  aperture  we  find  that  it  intercepts  an  aperture 
element 

^  =  - F -  di/  =  s,  di/  ,  (Cl) 

cos  i(/2 

where  is  the  distance  from  the  focus  to  the  reflection  point  on  the  cylinder. 
Thus,  the  ray  tube  has  the  same  cross-section  area  at  the  reflector  as  at  the 
aperture.  The  incident  field  E.(^)  at  the  reflector  is  assumed  normalized  such 
that  |E1I  is  the  intensity  in  Watts  per  meter.  Since  the  power  within  the  ray 
tube  is  constant  it  follows  that  the  normalized  aperture  distribution 


z 


where  the  phase  factor  expresses  the  phase  delay  over  the  distance  sa  from  the 
reflector  to  the  aperture. 

The  far  field  based  on  aperture  integration  is  now  obtained  as  follows.  Each 

aperture  element  dx  represents  a  source  of  strength  p(x)dx,  which  radiates  a 

(2) 

cylindrical  wave,  described  by  the  Hankel  function  H  (kr).  The  total  far  field 
C2  0 

can  be  shown  to  be 


Ea  =  f  p  (kr)  dx 


■-vS«'ikE  1  . 


where  the  superscript  a  denotes  aperture  integration,  and  we  have  made  the 

usual  far-field  approximations  for  the  Hankel  function  and  the  distance  r  from 

the  aperture  element  dx  to  the  observation  point.  In  the  sidelobe  region  the 

C3 

integral  reduces  essentially  to  two  end  contributions,  so  that 


-TT  -jkR  sinla-kDi 

~  -  Tli  p( d/2)  - .  _ 

jr  VlT  H  '  Sin© 


kD  sin©) 


Finally,  in  view  of  Eq.  (C2)  we  obtain 


a  a/n  sin(a-  kD  sin©) 

Ea  ~  -*Ji  £ _  E.(i/  )  - - 

ir  ^  i'v  o'  sm© 


The  far  field  E®  according  to  the  GTD  method  consist  of  the  sum  of  the  fields 
Ej  and  Eg,  scattered  at  the  left  and  right  edge  of  the  reflector,  as  shown  in 
Figure  C3.  Due  to  the  feed  symmetry  the  incident  field  is  equal  at  both  edges 
we  have  in  the  far  zone 


C2.  Harrington,  R.  (1961)  Time-Harmonic  Electromagnetic  Fields,  McGraw 
Hill,  N.  Y.  ,  p.  288. 

C3.  Felsen,  L. ,  and  Marcuvitz,  N.  (1973)  Radiation  and  Scattering  of  Waves, 
Prentice-Hall,  N.  J. ,  p.  387. 


where 


C  *  e”^*/4  sin(jr/n)  /  (nV2irk') 


n  =  2  -  t^/ir  . 


We  now  investigate  the  behavior  of  D1  +  D2  and  -  D2  in  the  sidelobe 
region  near  the  main  beam,  where  1 0 1 « 1.  In  general  n  and  it  will  be 
independent  parameters  and  therefore  we  obtain 


lim  (D.  +D„)  =  C 
0-0  1 


COS—  -  CO 

n 


s1q  sin2— 
n  n 


For  Dj  -  D2  however,  the  situation  is  different  and  there  we  find 


D.  -D,  ~  C  - - - — 

sin—  sin— 
n  n 


which  is  unbounded  when  0  —  0.  Substituting  the  above  expressions  into  Eq.  (C9) 
yields 


j*/ 


Isinu'  +  jc  cosu']  , 


where 


1  ,  cos  n  .  ir  .  6 

e  =  - —  + - = -  svn—  sin- 

cos— -  cos 2  sin  — 
n  n  n 


The  second  term  within  the  bracket  in  Eq.  (C13)  is  in  phase  quadrature  to 
the  first  one  and  has  a  much  smaller  magnitude.  Its  main  effect  is  to  fill  in  the 
pattern  nulls.  For  reflector  edge  angles  w^  <  90°  and  normal  F/D  ratios  we 


find  that 


|e(0,0o,n)|  <  I  sin— I 


Vv‘vv' 

■I-'vSV'1 

*  *  ?  a  .*r  v 


-  -N 


?v:v?jiv 

► 


’  ■  ■  ,  ’  .•.*.*«***ji**’i*.’i'*,»*»*t»*'»,»  ,J»1'  »  •  .  «  •  *  *  W>  I  «  «»»«*»»^ 


V  .  v 


1  ,>  «  » 

v'->» 

v\«V!<V 


The  second  term  may  therefore  be  neglected  and  as  long  as 


lei  s  e 


max  ”  4 


(C16) 


the  relative  error  in  the  field  pattern  envelope  will  be  less  than  1/32,  or  about 
3  percent.  Thus,  we  finally  obtain 


Eg-  ^sS/n)  eVr~W  sin(2  1 “  sine) 


Comparing  Eg  above  with  Ea  of  Eq.  (C 5),  we  find  that  the  ratio 

,2 


(C17) 


Si  =  nsin(^/n)  1+  A2  (l  1  \  <  !  +  =  1+i  . 

£g  slnlB  ®  \  n2  j  ®  96 


(C18) 


where  we  have  substituted  Eq.  (C16)  for  0 _ .  The  relative  difference  between 

max 

Ea  and  Eg  thus  is  l/96,  or  about  1  percent.  This  error  adds  to  the  3  percent 
relative  error  in  Ea  incurred  earlier. 

In  summary,  we  have  shown  that  the  sidelobe  pattern,  in  particular  the  pat¬ 
tern  envelope  and  the  sidelobe  positions,  is  accurately  determined  by  the  aperture 
integration  method  out  to  angles  of  ±  15°  or  ±  1/4  <D/X)  beamwidths  around  the 
main  beam.  In  this  sector  the  sidelobes  can  be  viewed  to  be  determined  by  the 
aperture  field  alone,  independent  of  the  reflector  edge  geometry.  Consequently, 
even  when  we  scan  the  main  beam  by  imposing  a  linear  phase  taper  over  the 
aperture,  we  can  expect  aperture  integration  to  accurately  predict  the  sidelobe 
pattern  over  a  ±l/4(D/X)  beamwidths  sector,  but  this  time  centered  around  the 
scan  direction. 

The  above  conclusions  are  independent  of  the  feed  polarization,  as  can  be 
easily  shown.  Also  they  apply  equally  well  for  a  circular  reflector  as  for  the 
cylindrical  reflector  considered  here,  since  the  diffraction  coefficients  are 
identical  in  the  two  cases.  However,  we  feel  that  the  idealized  geometry  con¬ 
sidered  here  may  be  a  particularly  favorable  case  and  that  caution  should  be 
exercised  when  applying  these  conclusions  to  more  general  problems,  as  for 
exa»  .pie,  offset  reflectors,  highly  tapered  illuminations,  or  illuminations  with 
large  amplitude  and  phase  errors. 
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